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We show that in second-order phase transformations in- 
duced by an inhomogeneous quench the density of topological 
defects is drastically suppressed as the velocity with which 
the quench propagates falls below a threshold velocity. This 
threshould is approximately given by the ratio of the healing 
length to relaxation time at freeze-out, that is at the instant 
when the critical slowing down results in a transition from 
the adiabatic to the impulse behavior of the order parameter. 
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In a uniform cosmological phase transition character- 
ized by an order parameter with a non-trivial homo- 
topy group, creation of topological defects is virtually 
inevitable. As pointed out by Kibble H some years ago, 
different domains shall select different low-temperature 
vacua, leading to irreconcilable differences which con- 
dense out into topological defects. Thus, special relativ- 
ity implies disorder in the context of Big Bang, yielding 
a useful lower bound on the initial density of defects. An 
analogous situation is also encountered in the condensed 
matter setting M, suggesting the possibility of experi- 
mental investigation of a cosmological scenario. For this 
case, however, a lower bound based on the light-cone 
causal independence is no longer useful. In condensed 
matter systems (and most likely also in the cosmological 
second-order phase transitions), the initial density of de- 
fects has to be computed through arguments which rely 
on the dynamics of the order parameter ||,|| ■ 

The estimate of defect density proposed by one of 
us relies on the observation that, as a result of critical 
slowing down, the state of the system which crosses the 
critical region at a finite pace will inevitably cease to 
keep up with the changes of thermodynamic parameters 
at some point sufficiently near the critical temperature 
In a homogeneous quench, this will happen every- 
where at the instant when the characteristic time e/e, at 
which changes in the mass term of the effective Landau- 
Ginzburg (LG) potential V{<j>) = (</> 4 - 2e(i)0 2 )/8 be- 
come faster than the relaxation time, i.e. r = To | e | — 1 
in the overdamped case. When e(t) ~ t/rQ, this leads 
to a freeze-out time t ~ (tqto) 1 ^ 2 . At t = — i the or- 
der parameter <j> becomes too sluggish to keep up with 
the evolving shape of the effective potential. An evo- 
lution dominated by dynamics driven by V(</>) shall re- 
sume t after the phase transition, which takes place when 
e(t — 0) = 0. The interval [— i, +i\ represents the im- 



pulse stage, namely the period during which the effect 
of dVdcj) is negligible, although fluctuations and damp- 
ing continue to matter. At t , e = e(t) = (to/tq) 1 / 2 , 
which leads to the estimate of the corresponding healing 
(or correlation) length £ = £ ( t q/ t o) 1 ^ 4 - This freeze- 
out healing length constrains the domain size relevant 
for defect formation. The corresponding relaxation time 
is f = (tqTo) 1 / 2 , and the characteristic phase ordering 
speed is v = £/f = v Q (tq/tq) 1 ^, where v = £gAb- 
This overdamped LG example can be generalized 

Homogeneous quenches are a convenient idealization 
and may be a good approximation in some cases. How- 
ever, in reality, the change of thermodynamic parameters 
is unlikely to be ideally uniform. Thus, the mass pa- 
rameter e(t, r), varying in both time and space, must be 
considered in defect formation. As a consequence, loca- 
tions entering the broken symmetry phase first could then 
communicate their choice of the new vacuum as the phase 
ordered region spreads in the wake of the phase transi- 
tion front. When inhomogeneity dominates, symmetry 
breaking in various, even distant, locations is no longer 
causally independent. The domain where the phase tran- 
sition occurred first may impose its choice on the rest of 
the volume, thus suppressing or even halting production 
of topological defects. 

Experiments carried out in 3 He 0B, where a small 
volume of superfluid is re-heated to normal state, and 
subsequently rapidly cools to the temperature of the 
surroundings, are a good example of an inhomogeneous 
quench: The normal region shrinks from the outside. 
Yet, topological defects are created, thus suggesting that 
the phases of distinct domains within the re-heated re- 
gion are selected independently. Even in 4 He, where the 
transition can be induced by a change of pressure, it is 
difficult to rule out the possibility that a quench may be 
somewhat non- uniform, thus causing decrease of the den- 
sity of defects, which could explain the recent evidence 
of non-appearance of vortices Q . 

Here we consider two idealized cases of inhomogeneous 
quenches: (i) a shock wave in which e(t, r) is a propagat- 
ing step in space and (ii) a linear front in which e(t, r) 
linearly interpolates between the pre- and post-transition 
values. For both scenarios, we investigate the threshold 
velocity v t at which phase ordered region expands be- 
hind the e = critical point. This velocity is given by 
the ratio of the healing length to relaxation time set by 
the dynamics. In an inhomogeneous quench, vt will be 
— as it was anticipated early on Q — of the order of v. 
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In an insightful paper directly stimulated by the 3 He 
experiments, Kibble and Volovik have argued that the 
initial density of defects should conform with the homo- 
geneous quench estimates of H,|| when the velocity v of 
the inhomogeneous quench front exceeds v; on the other 
hand, the initial density of defects should be suppressed 
as v/v for the case of a slowly spreading phase transition 
0. We find that v indeed plays a crucial role. How- 
ever, our study shows that the suppression is even more 
dramatic. When v > v, essentially no defects appear. 

SHOCK WAVE. To begin with, we consider a de- 
cay of a false symmetric vacuum to a true symmetry bro- 
ken phase in a one-dimensional dissipative 4 model. In 
dimcnsionlcss units, 



+ - 0" + - (0 3 -e o 0) = O 



(1) 



where 0(i, x) is a real order parameter and eo is a positive 
constant. We look for a solution (f>(t,x) which interpo- 
lates between 0(t, — oo) = —y/eo an d <j>(t, +oo) = 0. Such 
a solution can not be static; it is a stationary half-kink 



(f>{t,x) = -y/eS 1 + exp 



e (x - vft) 
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Our shock wave quench is a sharp pressure front prop- 
agating with velocity v. That is, 

4> + V<i> - 0" +~ (<f> 3 -e(t,x)<t>)=#(t,x) , (4) 



where 



e(t, x) = eo Sign(/j — x/v) 



(5) 



is the space-time dependent effective mass and §{t, x) is a 
Gaussian white noise of temperature 9 with correlations 



($(t,x) &(t',x')) = 2 ry 6 5(t - 1') 5{x - x') 



(6) 



There is a unique vacuum (0 — 0) to the right of the 
propagating front [x > vt), and the symmetry is broken 
(0 — ±y/eo) behind the front (x < vt). 

The field in the = vacuum does not respond in- 
stantaneously to the passing front. There are two quali- 
tatively different regimes: 

1) v > vf, the phase front propagates faster than the 
false vacuum can decay. The half-kink (|^) lags behind 
the front (^) ; a supercooled symmetric phase grows with 



velocity v 



The supercooled phase can not last for 



long; it is unstable, and the noise d makes it decay into 
one of the true vacua. Since the noise does not have any 
bias, the decay results in production of kinks. 



2) v < vf , the phase front is slow enough for a half- 
kink to move in step with the front, 0(i, x) — H v (x — vt). 
The symmetric vacuum decays into one definite non- 
symmetric vacuum, the choice is determined by the 
boundary condition at x — > —oo. No topological de- 
fects are produced in this regime. To make sure that the 
field can not flip, we must check if the stationary solution 
H v (x — vt) is stable against small perturbations. 

We investigate the stability in the r\ — * oo limit when 
the system is overdamped and the <j> term in Eq. (|]) can be 
neglected. H v is most likely to be unstable for v = v^~ . 
We use the ansatz 



rj v%x 



(f>(t, x) = H n (x - v s t t) + f(x - v s t t) exp Tt 



(7) 



The eigenequation turns out to be time-independent; 
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The "potential" in the square brackets is positive definite 
for = 3 v / eo/27/. This proves the stability of H v * at 
sufficiently low noise temperature 9. 

In the opposite 77 — > limit the half-kink (Q) just below 
the threshould uf « 1 becomes a step function H v * (x — 
Vft) s» y/eo[—l + Sign(x — vft)]/2. The potential on the 
r.h.s. of Eq.(j|) is again positive for any x. 

In summary, no topological defects are produced for 



v < vi 



At 



vf, there is a sudden jump in the 



density of defects left behind the shock wave. As v 
increases above w*, the density saturates at a value 
characteristic for an instantaneous uniform quench with 
e(t,x) — eo Sign(f). With increasing 9, the discontinu- 
ity at v = vf will be softened. For v ^> u|, where the 
quench is effectively homogeneous, the density of defects 
will grow logarithmically with 9 ]8j . 

These expectations are borne out by the numerical 
study of kink formation which uses the same code as in 
Rcf. |1. We illustrate them in Fig. [j] for 77 = eo = 1. 
For all but the highest temperature 9 = 0.1, there are 
essentially no kinks produced in quenches with the ve- 
locity of less than 0.8, which is in good agreement with 
our analytic estimate vf — 0.83. However, for the highest 
temperature, defects appear at a sub-threshold velocity. 
We note that at this temperature potential barrier sepa- 
rating the two minima of the LG potential is comparable 
with 9. 

LINEAR FRONT. Let us consider now the linear 
inhomogeneous quench 
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We assume that the linear part of e(t, x), namely the 
interval — sqvtq < x — vt < €qvtq, is much wider than the 
healing length to the left of the front, 26qvtq » I/^/cq. 
If not then the shock limit ([|) applies. 

In the absence of noise, the propagating linear front is 
followed by a stationary half-kink <p(t,x) = h v (x — vt). 
This half-kink lags a distance 5x behind the front. Sx 
can be estimated by similar arguments as those which 
led to Eq.(|3|). At Sx behind the front x — vt the mass 
parameter is eg x — Sx/vtq. The replacement eo —> £8x in 
Eq. (||) gives a velocity v t (Sx) the half-kink would propa- 
gate with if it were at Sx. The velocity increases with Sx. 
The half-kink gets stuck at such a Sx that this velocity 
is equal to the actual front velocity v, v t (Sx) — v. This 
takes place at Sx = 16t] 4 v 5 tq/81(1 — v 2 ) 2 , which grows 
quickly with v. 

When v is greater than u| in Eq.@, Sx > eqvtq and 
the half-kink does not stay in the linear regime. It enters 
the e = eo area where it moves forward with velocity 
wf < v. Like in the shock limit the supercooled phase 
grows at a constant rate and decays giving rise to kinks. 

When v < v", the half-kink remains confined in the 
linear regime. Even in this case, for v greater than certain 
threshould v t , the width Sx of the supercooled phase may 
be large enough for this phase to be unstable. A half-kink 
h v (x — vt) confined to the linear regime satisfies 



Defects can be produced for v > Vt, where 



(1-v 2 ) h'l{x)+r}vh' v {x) 



^-K(x)- 1 -hl(x) = o. 

(10) 



We rescale x — c\ x and h v = c 2 h in Eq.(|T(J), then set 
Ci = [2 v tq (1 — u 2 )] 1 / 3 and c 2 — y/ci/vTQ to obtain the 
rescaled equation 



L 1 [h](x) = h"+ph'-x h-h 3 = , 



(11) 



with primes now denoting derivatives w.r.t. x. Eq. (|1 
has a single parameter 



P 



2 l/3 ??u 4/3 T l/3 
(1 - V 2 ) 2 / 3 



(12) 



The half- kink h becomes unstable at a threshold p = p t . 
At this critical pt, h has a zero mode Sh, which satisfies 



L 2 [Sh](x)=Sh +p t Sh -xSh-ih 2 Sh = . (13) 

The value of pt was found in two steps. First, we found 
solutions to Eq. ([n]) for a range of p with the relax- 
ation scheme h = Li[h]. We applied then the relaxation 
scheme Sh — L,2[Sh] to the linear Eq. ( |l3| ) with the ini- 
tial condition 8h(t — 0,x) = 1. The field relaxed to 
Sh(t — > oo,x) = for p < p t , and it blew up without 
limit for p > p t . For p « p t , we observed a long lived 
localized zero mode structure. The threshold estimated 
in this way is pt = 6.5 . . . 6.6 . 
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The instability appears because the eigenvalue of the 
lowest mode of a linearized fluctuation operator around 
h v passes through zero when v = v t . The passage is 
smooth, so we do not expect any discontinuity in the 
density of defects as a function of v. For the same reason, 
we expect the threshold at v t to be gradually softened 
with increasing noise temperature 9. For i; » «[, the 
inhomogeneity of the quench is irrelevant, and the density 
of defects can be estimated by scaling argument HQ as 
for a homogeneous quench with a timescale tq. 

This analysis is confirmed by the numerical study of 
linear quenches shown in Fig. ^|. For the lowest temper- 
atures, there are no kinks formed below the threshold, 
which for our 77 = 1 is vt ~ 0.77. However, as tempera- 
ture increases from 9 = 10~ 10 to 9 = 0.1, kinks begin to 
appear at velocities as low as ^0.42. This decrease of the 
threshold for kink formation is now more gradual than 
for the shock wave case of Fig. [j]. 

CONCLUDING REMARKS. We found that for 
both the shock and the linear inhomogeneous quench, 
the density of defects is drastically suppressed for quench 
velocities lower than the characteristic velocity Vt ~ v. 
This prediction was verified by numerical simulations for 
kinks in one dimension. The theory can be generalized 
to higher dimensions and to complex order parameter in 
a straightforward manner. 

More quantitative understanding of the dependence of 
the number of defects on the quench velocity requires 
further investigation. The non-hermitian nature of the 
operator L 2 [^] in Eq.(|l3|) should be carefully taken into 
account. 

Our prediction that no defects are produced below cer- 
tain threshold is in contrast with Ref. |t]], where some 
defects are predicted even at very low velocities (but see 
Ref. [[[o| for a different conclusion). Formation of planar 
solitons in 3 He — A JOL which takes place at velocities 
10 2 — 10 3 times lower than the threshold velocity, seems 
to confirm the latter theory. Let us however make two 
remarks: 

1) No uniform temperature gradient was deliber- 
ately applied in this experiment although some non- 
uniformities could result in local temperature gradients. 
There may have existed such regions of the sample where 
the gradient was null or negligible. Solitons could be eas- 
ily created in these areas. This explanation can be ver- 
ified in the same experimental set up; application of a 
sufficiently strong uniform temperature gradient across 
the sample should suppress any soliton formation. 

2) Moreover, in three dimensions, the situation is 
slightly more complicated. Straightforward generaliza- 
tion of our analysis would rule out formation of zero- 
dimensional defects, as well as one and especially two 
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dimensional defects parallel or askew to the direction of 
the quench front. However, such extended string and 
domain-wall like structures could still "grow" in the di- 
rection perpendicular to the front (and parallel to the 
direction of its propagation), providing that their seeds 
exist, say, along the wall where the symmetry is broken 
first. Strings would grow from some seeds and antistrings 
from the other. The growth of individual (anti-)strings 
would not be perfectly perpendicular to the front, they 
would be wandering around in a chaotic manner. From 
time to time the growing ends of a string and an anti- 
string would meet and coalesce to form a "jumping rope" 
with its other ends anchored at the original seeds. For 
global strings (like vortices in 4 f/e) the coalescences are 
accelerated by a long range mutual string-antistring at- 
traction. Such a rope would shrink to the original wall 
thus disappearing from the bulk. At some stage all possi- 
ble coalescences would have already taken place leaving 
only a net number of say, strings, proportional to the 
square root of the number of seeds. Some of these sur- 
vivors would be forced by their mutual repulsion to ter- 
minate on the side walls. Only a fraction of them would 
reach the opposite wall spanning through the bulk where 
they can be unambiguously detected. Their density is 
likely to be, in any case, orders of magnitude below the 
estimates based on [^]||. Similar argument applies to 
membrane-like solitons in |Tl| ] . 

We think that inhomogeneity is a factor which may 
also need to be taken into account in interpretation of 
the recent negative 4 He experiment ||. 
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FIG. 1. Density of kinks n as a function of velocity v for 
the shock wave (0) with 77 = eo = 1 (overdamped system). 
In this overdamped regime, the predicted threshold velocity 
is vt ~ 0.83. The plots from top to bottom correspond to 
6 = 10~\ 10" 2 , 10" 4 , 10" 6 , 10~ 8 , 10~ 10 . At low 6, we get a 
clear cut-off velocity at v « 0.8, which is consistent with the 
prediction. 



[1] T.W.B. Kibble, J. Phys. A9, 1387 (1976); Phys. Rep. 

67, 183 (1980). 
[2] W.H. Zurek, Nature 317, 505 (1985). 
[3] W.H. Zurek, Acta Phys. Polon. B24, 1301 (1993); Phys. 

Rep. 276, 177 (1997). 
[4] V.M.H. Ruutu, V.B. Eltsov, A.J. Gill, T.W.B. Kibble, 

M. Krusius, Y.G. Makhlin, B. Placais, G.E. Volovik and 

Wen Xu, Nature 382, 334 (1996); 
[5] C.Bauerle, Yu.M. Bunkov, S.N. Fisher, H. Godfrin and 

G.R. Pickett, Nature 382, 332 (1996); 
[6] M.E. Dodd, P.C. Hendry, N.S. Lawson, P.V.E . McClin- 
tock and C.D.H. Williams, |cond-mat/9808117 . 
[7] T.W.B. Kibble and G.E. Volovik, Pis'ma v ZhETF 65, 

96 (1997). 



4 



0.5 



V 



1 



1.5 



FIG. 2. Density of kinks n as a function of velocity v 
for the linear inhomogeneous quench Q with tq — 64 and 
r\ = 1. The predicted threshold is v t = 0.77. This cut-off is 
achieved for low 0. The plots from top to bottom correspond 
to 6 = 10" 1 , 10" 2 , 10" 4 , 10" 6 , 10" 8 , 10" 10 . 
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